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ABSTRACT 



Recent surveys have identified seven hypervelocity stars (HVSs) in the halo of 
the Milky Way. Most of these stars may have originated from the breakup of binary 
star systems by the nuclear black hole SgrA*. In some instances, the breakup of the 
binary may lead to a collision between its member stars. We examine the dynamical 
properties of these collisions by simulating thousands of different binary orbits around 
SgrA* with a direct N-body integration code. For some orbital parameters, the two 
stars collide with an impact velocity lower than their escape velocity and may therefore 
coalesce. It is possible for a coalescing binary to have sufficient velocity to escape the 
galaxy. Furthermore, some of the massive S-stars near Sgr A* might be the merger 
remnants of binary systems, however this production method can not account for most 
of the S-stars. 

Key words: black hole physics-Galaxy:center-Galaxy:kinematics and dynamics- 
stellar dynamics 



1 INTRODUCTION 



First theorized bv lHillsl (llOSSl ). a hypervelocity star (HVS) 
has sufficient velocity to escape the gravitational pull of the 
Milky Way galaxy. The first HVS, SDSS J0907 45.0-H024507, 
was recently discovered i n the Galactic halo (iBrown et al.l 
l2005l : iFuentes et al.l [20051 '). This HVS is located at a helio- 
centric distance of ~ 110 kpc and has radial velocity 853 
± 12 kms~^, over twice that needed to escape the gravita- 
tional pull of the Milky Way. Since th at initial discovery , 
six other HVSs have been identified ( l Edelmann et al.|[2005l : 
Hirsch et aLll2005l : iBrown et ai1l2006al : iBrown et al.ll2006bh . 



Hillsl (|l988l ) suggested that a HVS might result from a close 



encounter between a tightly bound binary star system and 
the bl ack hole at the Galactic center, SgrA*. lYu fc Tremaind 
l|2003ll refined Hills' argument and noted that HVSs might 
also be produced by three-body interactions between a star 
and a binary black hole system. Because the existence of a 
second (intermediate-mass) bla ck hole in the Galactic center 
(iHansen fc Milosavlievicll2003l ) is only a hypothetical possi- 
bilitv jSchodel et aLbooa T we focus here on the disruption 



of a tightly bound binary by a single supermassive black 
hole (SMBH) with a mass of ~ 4 x 10® M© whose exis- 
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tence is robustly supported by data (e.g. Ghez et al] l2005l : 
iReid fc Brunthaleij|2004l : ISchodel et alll2003l ). 



Simulations show that tight binaries can produce HVSs 
with velocities compa rable to the observe d HVSs (e.g. 



Ginsburg fc Loebl |2006| (hereafter Paper I) ; iBromlev et al.l 



20061). In Paper I, we show that the companion to the hy- 



pervelocity star will be left in a highly eccentric orbit, which 
agrees with the known orbits of a number of S-stars orbit- 
ing Sgr A* (e.g. Ec kart fc Genzelll 19971 . ISchodel et al]|2003l . 
and lGhez et al...2005l )7 Therefore, we suggested that some of 
these stars are former companions of HVSs. Furthermore, a 
small fraction (~ 10%) of the binary systems were found to 
collide. Here we examine in detail the dynamical properties 
of such collisions, and check whether some of these collisions 
may end in coalescence. 



In §2 we describe the N-body code and the simulation 
parameters that were used. In §3 we discuss our numerical 
results for the collisions, and in §4 we discuss the outcome 
of binary mergers at the Galactic Centre. Our goal is not to 
cover the entire range of binaries that could produce HVSs 
or end in collisions, but rather to determine whether some 
tight binaries with masses similar to the HVSs observed thus 
far could coalesce. 
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2 COMPUTATIONAL METHOD 

In our s tudy we have used the N-body code written by 
Aarseth (|Aarsethlll999l ') whose details were described in Pa- 
per I. We treat the stars as point particles and ignore tidal 
and general relativistic efltects on their orbits, since these ef- 
fects are small at the distance ('^ lOAU) where the binary 
is tidally disrupted by the SMBH. We have set the mass of 
the SMBH to M = 4 x W^Mq. The masses of the binary 
members are set to either 3Mq & 3M© [sin ce SMq is the 
estim ated mass of SDSS J090745. 0-1-024507 (iFuentes et alj 
120051 )] ■ or to 3Mq & lOM© [as l OMa is comparable to t he es- 
timated mass of HE 0437-5439 (|Edelmann et al.ll2005l ')]. Ah 
runs start with the center of the circular binary located 2000 
AU (= 10~^pc) away from the SMBH along the positive y- 
axis. This distance is comparable to the inn er scale of the ob- 
served distribution of stars around SgrA* (lEckart &: Genzell 
ll997l : ISchodel et al.ll2003l : iGhez et al.ll2005^ . allowing the re- 
maining star to populate this region after the ejection of its 
companion. This radius is also much larger than the binary 
size or the distance of closest approach necessary to obtain 
the relevant ejection velocity of HVSs, making the simulated 
orbits nearly parabolic. 

We ran two sets of data. The first had the binary system 
rotating along the x~y plane and the second along the y- 
z plane. We used the same initial distance for all runs to 
make the comparison among them easier to interpret as we 
varied the distance of closest approach to the SMBH or the 
relative positions of the two stars within the binary. We 
chose initial binary separations between a — 0.05 and 0.2AU 
because such a range is likely to produce HVSs for the above 
parameters (see Paper I). Sig nificantly wider binarie s would 
give lower ejection velocities (|Gualandris et al]l2005h . Much 
tighter binaries would not be easily disrupted by the black 
hole, or may coalesce to make a single star before interacting 
with the SMBH. The radius of a main sequence star of a few 
solar masses is ~ O.OIAU, and that of a 10 solar mass star is 
~ 0.03AU (see, e.g. Fig. 4 in Freitag & Benz 2005). Binaries 
tighter than ~ 0.02AU are precluded because the two stars 
will develop a common envelope and eventually coalesce. 

In the Galactic disk, about one-third to half of all stars 
form in binaries or small multip le systems (see e.g. iLadal 
I2OO6I : [Puauennov fc Mavorlll99ll ). with roughly equal prob- 
ability per logarithmic interval of separations, dP/dln{a) = 
const (e.g. ^t 1983; Hcacox 1998 ; Larson 2003). In the 
Galactic center environment, the maximum binary separa- 
tion is limited by the tidal force of SgrA* at the distance d 
where the binary is formed (fo r conditions that enable sta r 
formation near the SMBH, see iMilosavlievic" fc Loeblliooi ). 
Since the mass of the black hole is ~ 10® times larger than 
that of a star, this implies a maximum binary separation 
less than (10"'')^''^ = 10"^ of the initial distance d. For 
d = 2 X 10"^ AU, the upper limit on the binary separation 
would be 20AU (or smaller if the tidal restriction applies 
during the formation process of the binary). If we assume 
a constant probability per ln(a) for 0.02 < a < 20AU, then 
the probability of finding a binary in the range of a = 0.05- 
0.2AU is substantial, ~ 20%. 

As shown in Paper I, the initial phase of the binary orbit 
plays a crucial role in the outcome. Therefore, we sampled 
cases with initial phase values of 0-360 degrees in increments 
of 15°. As initial conditions, we gave the binary system no 



radial velocity but a tangential velocity with an amplitude 
in the range between 5 and 25 km s~^ at the distance of 
2000AU. In total, we ran 2000 simulations. 



3 PROPERTIES OF HYPERVELOCITY 
COLLISIONS 

Given a binary system with stars of equal mass m separated 
by a distance a and a SMBH of mass M S> m at a distance b 
from the binary, tidal disruption would occur if & < 6t where 



M 



(1) 



The distance of closest approach in the initial plunge of the 
binary towards the SMBH can be obtained by angular mo- 
mentum conservation from its initial transverse speed v± at 
its initial distance from the SMBH, d, 



V_Ld - 



1/2 



(2) 



The binary will be tidally disrupted if its initial transverse 
speed is lower than some critical value. 



1/2 



10^ 



,1/6 



ds.a 



km s 



(a/0.1 AU), da.s = (d/2000 AU), mo.: 



(3) 



where a_] 

(m/3M0). If v± < u_L,crit, one of the stars receives sufficient 
kinetic energy to become unbound, while the second star is 
kicked into a tighter orbit around the SMBH. The ejection 
speed Vcj of the unbound star can be obtained by considering 
the change in its kinetic energy ~ v5v as it acquires a ve- 
locity shift of order the binary orbital speed Sv ~ ^ Gm/a 
during the disruption process of the binary at a distance 
~ bt from the S MBH when the binary cente r-of-mass speed 
is w ~ ^GM/bt (|Hillsll 19881 : IYu fc Tremaing,2003 ). At later 
times, the binary stars separate and move independently 
relative to the SMBH, each with its own orbital energy. For 
V < wx,crit, we therefore expect 



(Gm\' (GM\ 



1/2 



1/2 



= 1.7 X lO'^m, 



3 1/3 



'1/2 
-1 



km s 



(4) 



Under some circumstances, the binary is disrupted in 
such a way that the two stars collide. Assuming that the 
impulsive kick is given by the SMBH towards a random di- 
rection within the orbital plane and ignoring gravitational 
focusing (which is important at low speeds), the probability 
for a collision in a case that otherwise would have produced 
a HVS is four time the radius of a star divided by the cir- 
cumference of a circle with a radius equal to the binary 
separation. The likelihood for a collision is expected to be 
smaller in the more general case where the binary lies in a 
different plane than its orbit around the SMBH, unless grav- 
itational focusing dominates. Table [T] summarizes the actual 
statistical results from our runs. 

The two stars would merge as a result of the collision 
if their relative speed is lower than the escape speed from 
their surface (~ 500 km s~^). In our runs 22% of all col- 
lisions have impact velocities low enough to allow the two 
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Figure 1. Fraction of all collisions (in percent per 100 kms~^ bin) versus relative velocity upon impact (in km s~^). The left section 
is for amin = 0.02 AU and the right section is for a^in = 0.04 AU. The label of the lower left panel corresponds to all panels. The 
dashed vertical line shows the impact velocity that would have resulted from free fall starting at the binary separation (see equation 
(GJ. The solid line is the median velocity of all runs. (We choose to use the median rather than the average value because outliers bias 
the data otherwise.) The minimum impact parameter for a collision is expected to be amin = {Ri + R2) = 0.02 AU for a SMq &c 3Mq 
binary, but amin = 0.04 AU for a IOMq & SMq binary. We show results in other cases for pedagogical purposes, namely to illustrate 
the dependence of the results on the binary masses and a^in separately. If the IOM0 companion is a black hole then amin ~ 0.01 AU. 



a (AU) 


P(3Mq) 


P(1OM0) 


0.05 


0.11±0.02 


0.21±0.05 


0.10 


0.11±0.02 


0.13±0.04 


0.15 


0.06±0.01 


0.12±0.03 


0.20 


0.03±0.01 


0.04±0.02 


0.05 


0.09 


0.27 


0.10 


0.06 


0.13 


0.15 


0.03 


0.09 


0.20 


0.02 


0.07 



Table 1. Collision probability with different values of a for bina- 
ries of 3M0 & 3M0 (second column) and lOM© & 3M0 (third col- 
umn). The top four rows show the values obtained from our sim- 
ulations with their corresponding Poisson errors. For comparison, 
the bottom rows show the expected probability from a simplis- 
tic "billiard ball" model (without gravitational focusing) in which 
the probability of a collision is 2{Ri + R2) /2TTa. Here {Ri}i— 1^2 
are the radii of the two stars and a the binary separation. 



stars to coalesce (see Table [2]). Also of note is the fact that 
many collisions involve hypervelocities of w > 1000 kms"'^ 
upon impact. The typical impact velocity of the two stars 
can be crudely estimated from a model in which the SMBH 
removes the angular momentum from the binary and causes 
the two stars to fall toward each other from their initial 



orbital separation. Conservation of energy 

Gmim2 



mim2 .2 
■r 



2 mi + m2 
yields the relative velocity upon impact 

1 



= const. 



Vf = 



2G{mi +m2){- 



-.1/2 



(5) 



(6) 



The actual impact speed would vary around this value due 
to the additional velocity induced by the SMBH tidal force 
along the axis connecting the stars. Nevertheless, Equation 
[6] agrees well with the median of the distribution of impact 
speeds in our runs (see Figure [1]). Collisions always occur 
shortly after tidal disruption, as seen from the separation of 
the black and filled circles in Figure (2] 



4 FATE OF THE COALESCING BINARY 

Stellar collisions ar e likely the ma i n assembly line of blue 



stragglers (see e.g. LeonardI Il989l : iBailvn fc Pinsonneaultl 



ll995l : lLombardi et al.ll2002l). and ultracompact X- ray bina- 
ries (e.g. Ivanova et al. 20051 : iLombardi et al.|[2006l ) in glob- 



ular clusters. The Galactic Centre of the Milky Way is an- 
other place where collisions are likely to occur. Tidal disrup- 
tions of a binary by the SMBH will produce ~ 0.1 collisions 
per HVS (see Paper I). The ultimate fate of the binary de- 
pends on the velocity of its member stars upon impact. As 
evident from Figure [1] the impact velocity Uimp can vary 
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Figure 2. Orbits of stars (in units of AU) prior to collision for binaries of SMq & 3Mq in panels (a) — (c) and binaries of IOMq & SMq 
in panels (d) — (/). For all panels, the graphs on the left are the orbits of the stars as they pass near the SMBH located at the origin. 
For panels (a) — (c), the graphs on the right are plotted at the center of mass frame. For panels (d) — (/), the graphs on the right are 
plotted at the rest frame of the IOM0 star. The red filled circle denotes the collision instant and the black filled circle denotes the time 
when the stars started to move towards each other as a result of the SMBH tidal force. Note that before the binary is disrupted, the 
3M0 star makes many revolutions as denoted by the blue circles. After the binary is disrupted, the approaching stars are denoted by 
black triangles for clarity. Panels a and d both have an initial separation of a = 0.05 AU. Panels b and e have a = 0.10 AU, and panels 
c and / have a = 0.15 AU. 



a (AU) 


P(3M0) 


P(1OM0) 


0.05 


0.46 


0.13 


0.10 


0.19 


0.05 


0.15 


0.08 


0.00 


0.20 


0.00 


0.00 


0.05 


N/A 


N/A 


0.10 


0.81 


0.33 


0.15 


0.93 


0.05 


0.20 


0.25 


0.03 



Table 2. Probability of coalescence upon collision for different 
values of a for binaries of 3M0 & 3M0 (second column) and 
IOM0 &i 3M0 (third column). The top four rows show the val- 
ues obtained from our simulations for amin = 0.02 AU, and the 
bottom rows for ami„ = 0.04 AU. Here we assume that the two 
stars will merge if Vimp < 500 km s~^. 



over a wide range of values. A star with Vimp > Vesc = 
[2G(mi + m2)/(7?i + R2)Y^^ will likely p ass through the 
other star, even during a head-on collision l|Freitag fc Bend 
However, in any collision there certainly will be in- 
teractions where the smaller st ar may gain mass an d the 
larger star will likely lose mass (|Freitag fc Benj|2005l ). Fur- 
thermore, a collision where the impact velocity is less than 
the escape velocity v^sc will not necessarily end as a merger. 
A grazing collision might result in envelope-ejection but no 



core merger, whereas a head-on collision mig ht result in core 
merg er, and thus form a more massive star ijPale fc DaviesI 
I2OO6I ) . The results of smoot hed particle hydrodynamics sim- 
ulations of blue stragglers (|Sills et al.l [200S) show that off- 
axis collisions initially have large angular momentum but 
eventually lose it to allow the merger to contract down to 
the main sequence. OfT-axis collisions, which are more proba- 
ble than head-on collisions, could nevertheless lead to HVSs 
which are rapidly spinning ( Alexander fc Kumar 2001). Fi- 
nally, a merger between a lower mass star with a higher mass 
star may extend the massive star's main-sequence phase 
l|Dale fc Daviesll2006h . 

lEdelmann et al.1 (|2005l ) notes that HVS HE 0437-5439 
might be the merger of two 4M0 stars, and that such a 
merger is consisten t with the age of the HVS. Furthermore, 
iHirsch et all (|2005l ) suggests that HVS US 708, a sublumi- 
nous O star, might be the merger of two helium-core white 
dwarfs. After losing mass, some of the coalescing binaries 
in our runs might end up with sufficient velocity to escape 
the galaxy. Unfortunately, our simulations cannot treat mass 
loss as well as the lar ge amount of thermal e nergy deposited 
in each collision (see [Leonard fc Livioll 19951 ). Coalescing bi- 
nary systems that remain bound to the SMBH could end up 
as massive S-stars. 

Aside from stellar collisions of main-sequence stars, it 
is possible for collisions to involve compact objects. As long 
as the colliding objects are gravitationally bound, the com- 
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pact object will eventually settle to the center of the merger 
remnant due to angular momentum transport by dynami- 
cal friction (gravitationally-induced spiral arms) or viscos- 
ity in the stellar envelope ( see e.g. the simulation of a black 
hole-helium star merger in IZhang fc Fr"TOill200ll) as well as 
gravitational radiation. The situation of a stellar mass black 
hole surrounded by a star also appears in the coUapsar 
prog enitors of long-duration Gamma-Ray Bursts (GRBs) 
(|MacFadven fc Woosley , .1999il . However, these events result 
from the collapse of the stellar core and so the accretion rate 
into the black hole is larger than the Eddington limit by 
more than 10 orders of magnitude (jNaravan et al.|[200ll ). In 
our case, the accretion might be limited by the Eddington lu- 
minosity and so the resulting sources would be much fainter 
than GRBs. It is unclear whether this accretion would result 
in an implosion or an explosion. In the case of a neutron star 
companion, on e gets a Thorne-Zytk ow object, with a similar 
accretion rate (|Podsiadlowskilll996l l. 



5 CONCLUSIONS 

Our N-body simulations show that tight binaries with a sep- 
aration between 0.05-0.20 AU which approach within a dis- 
tance ~ 10 AU from SgrA* could produce both HVSs and 
collisions (in ~ 10% of all cases) where occasionally the two 
objects may coalesce. Coalescence occurs if th e impact veloc- 
ity is sufficiently low (|Freitag fc Benj|2005l ). otherwise the 
SMBH would likely eject the two stars into separate orbits 
(see Paper I). The large variance in collision velocity shown 
in Fig. [T] will result in a broad range of collision products. 

Given the low production rate of HVSs in the Milky- 
Way galaxy ~ 10"^ yr"^ (Brown et al. 2006b), and that 
collisions occur in ~ 10% of all cases with mergers account- 
ing for ~ 20% of all collisions, the expected production rate 
of S-stars from binary mergers is once every 5 million years. 
Since this period is comparable to the lifetime of the massive 
S-stars, most of these stars can not be the merger remnants 
of binary systems. It is unlikely that any of these collision 
events will be observed in real time. Following a coalescence 
event, the remnant will appear different from main-sequence 
stars of the same mass only if the progenitors evolved sub- 
stantially before their collision time. 
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